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This study presents a systematic and quantitative analy- 
sis of the effect of inhomogeneous surface albedo on short- 
wave cloud absorption estimates. We use 3D radiative trans- 
fer modeling with gradually complex clouds over a simpli- 
fied surface to calculate cloud absorption. We find that 
averaging surface albedo always underestimates cloud ab- 
sorption, and thus accounting for surface heterogeneity al- 
ways enhances cloud absorption. However, the impact on 
cloud absorption estimates is not enough to explain the dis- 
crepancy between measured and model-calculated shortwave 
cloud absorptions. 


1. Introduction 

The anomalous shortwave cloud absorption, defined as 
the difference between measured and model-calculated ab- 
sorptions, has been found on the order of 10 W/m 2 . Re- 
gardless of the recent intense debates about whether this 
difference is significant [Ackerman et ah, 2003; O’Hirok and 
Gautier , 2003; Valero et a/., 2003], there is no doubt that 
some degree of discrepancy exists in observed and calcu- 
lated cloud absorptions, and this difference tends to be a 
bias rather than a random error [Valero et ah, 2003]. This 
excess absorption is not large but any such bias is of concern 
since radiative transfer models are widely assumed to be suf- 
ficiently accurate for climate modeling and remote sensing 
applications. 

Numerous efforts have been made to identify potential 
sources of this shortwave cloud absorption bias, including 
from the aspects of sampling issues, measurements uncer- 
tainty, cloud inhomogeneity, microphysics optical proper- 
ties, and aerosols loading [ Barker , 1992; Marshak et a/., 
1997; Valero et a/., 1997a; Cess et ah, 1999; Knyazihhin 
et ah, 2002; Ackerman et ah, 2003; O’Hirok and Gautier, 
2003; Oreopoulos et ah, 2003]. Based on high-resolution 
spectral albedo data along with a state-of-the-art radiative 
transfer model [Li et a/., 2002], it also has been stated [Li 
et ah, 2003] that accounting for the heterogeneity of sur- 
face albedo could reduce the systematic difference between 
measured and modeled cloud absorptions. However, the in- 
fluence of inhomogeneous surface albedo, if any, has been 
ignored in most climate models. As a result, up to now, 
there have been no systematic and quantitative analyses of 
the effects of surface heterogeneity on cloud absorption es- 
timates. This study aims to provide a thorough analysis 


to understand how accounting for surface heterogeneity af- 
fects cloud absorption, and to examine whether the observed 
cloud absorption bias could be explained by inhomogeneous 
surface albedo. 

2. Approach 

We used the Discrete-Ordinate-method radiative trans- 
fer model (DISORT) [Stamnes et ah, 1988], a Monte Carlo 
method [Marchuk et al ., 1980], and the Spherical Harmonics 
Discrete Ordinate Method (SHDOM) [Evans, 1998] to calcu- 
late cloud absorption. Models were set up with clouds over 
a checkerboard albedo surface (shown in Fig. 1), where the 
complexity of clouds increases from homogeneous to bro- 
ken clouds. The checkerboard surface was changed from 
the extreme case of black and white albedos, having the 
largest contrast, to a black and gray pattern, which is closer 
to measured albedos from the Atmospheric Radiation Mea- 
surement (ARM) Program. Cloud properties are defined via 
cloud optical depth r and single scattering albedo zoo, and 
cosine of the solar zenith angle (SZA) is denoted as fiQ. 



Figure 1. Schematic illustration of model setup. 

Based on energy conservation, cloud absorptance A can 
be computed from reflectance R and transmittance T as 

A(a) = 1 - R(a) - (1 - ct)T(a), (1) 

where A , R , and T are all functions of Lambertian surface 
albedo a, and (1 — a)T presents total surface absorption. 
Note that A, R, and T are also functions of r, zoo and po- 
Since [Liou, 2002] 
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two ends is around 8%. When the scale ratio decreases to- 
wards zero, the cloud absorptance is approaching the SIPA 
estimate. This limiting case of small s, in which the cloud 
layer is low, and thus most cloud absorption is attributed to 
photons directly from the underneath surface pixel, is close 
to the surface independent pixel assumption. On the con- 
trary, cloud absorptance reaches the homogeneous-surface 
estimate in the limiting case of larger s. Since this case cor- 
responds to a higher cloud base, photons from other surface 
pixels would have more chances to enter the cloud pixel of 
interest. As a result, the surface looks gray to the cloud, 
which is close to the homogeneous-surface assumption. 


4- Inhomogeneous clouds 

The previous section has demonstrated the effects of sur- 
face heterogeneity on cloud absorption for homogeneous 


Surface: Black and white checkerboard 
Cfoud: Homogeneous 
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Figure 4. Cloud absorptance as a function of scale ratio 
s for the ”biggest~effect” case: (r, za 0i po) = (16, 0.99, 
1.0). s is the ratio of cloud base height to the scale of in- 
homogeneous albedo surface. SIPA indicates that cloud 
absorptance approaches the estimate made by the sur- 
face independent pixel assumption in the limiting case of 
small $. HS shows that the cloud absorptance is toward 
the homogeneous-surface estimate in the limiting case of 
large s . 
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Figure 5. Cloud absorption as a function of scale ratio 
with an underlying black and white checkerboard surface 
for (a) the case same as Fig. 4, and (b) the case integrat- 
ing over various clouds, wavelengths, and solar zenith 
angles. The percentage indicates the relative bias in per- 
centage between the smallest and largest scale ratios. 


clouds. However, observed clouds are rarely uniform without 
any internal or horizontal variability [ Davies , 1978]. Hence, 
we need to understand whether inhomogeneous albedo sur- 
face affects heterogeneous clouds in the same way and with 
a comparable degree of impact. To simulate cloud inhomo- 
geneity, a fractionally integrated cascade model [Schertzer 
and Lovejoy , 1987] was used to generate various cloud struc- 
tures. For a variety of cloud fractions 50 realizations were 
produced at a given mean optical depth of 16. These re- 
alizations then were used to calculate cloud absorption for 
zdq = 0.99 and various zenith angles and cloud scales. We 
find that there is no significantly qualitative or quantitative 
difference in the effects of surface heterogeneity for fractal 
clouds rather than homogeneous clouds (figures not shown) . 
In 3D clouds, 5^7% of change in cloud absorption is found 
between the two limiting cases of SIPA and HS. 


5. Inclusions of broadband spectrum, diurnal 
cycle and various clouds 


We have demonstrated how inhomogeneous surface af- 
fects single- wavelength cloud absorption for homogeneous 
and heterogeneous clouds over a black and white checker- 
board surface. When we consider all variability of r, zu o, 
and fio , by defining cloud absorptance as A(r, tz?o, A^o, c*)> 
the overall cloud absorption (A) can be computed by 


{A) = j / A{t,w 0 ,Mo,a)- 

J ts?o * MO ^ T 

Pr(T)P, 10 (m))Pv,o(‘ cu o)dTdfiodwo, (7) 


where Fr , Pp.Q y and P«, 0 represent the probability density 
functions of r, po, and U7o> respectively. P T was estimated 
from simulations of the cascade model with a sample mean of 
16. P M0 was approximated by equal weights at zenith angles 
of 30, 45 and 60°. For Pt» 0 , we divided the solar spectrum 
into five intervals to find weighting factors for the corre- 
sponding averaged single scattering albedo. Then, using a 
simple quadrature rule, the approximated cloud absorption 
could be obtained. 

The resulting broadband cloud absorptance over a black 
and white checkerboard surface is shown in Fig. 5(b). For 
comparison, the biggest-effect case for homogeneous clouds 
is plotted in Fig. 5(a). The bias due to the change of 
scale ratio reduces to less than 4% after considering differ- 
ent clouds, a diurnal cycle, and a broadband spectrum. Note 
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Figure 6. Same as Fig. 5, but with a black (albedo of 0) 
and gray (albedo of 0.5) checkerboard surface in which 
the albedo is closer to the ARM measurements. 
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Popular Summary 


“Enhanced shortwave cloud absorption” (the discrepancy between measured and model- 
calculated cloud absorptions) has been a major concern in the climate community. The 
reason is that this excess absorption is always a bias on the order of 10 W/m 2 rather than 
a random error, and in turn, may have significant impact on climate modeling and remote 
sensing application. Among various explanations for this bias, it has been stated that 
taking surface heterogeneity into account could eliminate this discrepancy. However, the 
influence of inhomogeneous surface albedo, if any, has been ignored in most climate 
models. As a result, up to now, there have been no thorough analyses of the effects of 
surface heterogeneity on cloud absorption estimates. 

This study presents a systematic and quantitative analysis of the effect of 
inhomogeneous surface albedo on shortwave cloud absorption estimates. We provide 
both theoretical proof and numerical calculations to demonstrate that the use of an 
averaged surface albedo always underestimates cloud absorption. Thus, accounting for 
surface heterogeneity always enhances cloud absorption. However, the impact of 
inhomogeneous surface albedo (less than 0.5%, or ~1 W/m 2 ) is negligible and cannot 
explain the enhanced shortwave cloud absoiption. 



